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Abstract 

A strong magnetic field can make it advantageous to work in a coordinate system aligned with dipolar 
field lines. This monograph collect the formulas for some of the most frequently used expressions and 
operations in dipole coordinates. 

In some physical systems (e.g., the terrestrial ionosphere or the solar corona) the magnetic field can 
impose a strong anisotropy by restricting transport processes perpendicular to the lines of force. If the field 
is approximately dipolar it can be useful to work in an aligned coordinate system even though the vector 
operations are somewhat more complicated than in a Cartesian or spherical polar representation. Here 
we present several of the formulas that frequently occur when working with vectors in dipolar coordinates. 
Although some of the results have previously appeared in the literature [21 E] , others seem to be new. 

There are many possible (related) choices for the dipolar coordinates. The right-handed orthogonal 
system considered here, (q,p, <j)), is defined in terms of the usual spherical polar coordinates, (r, 9, 4>), by 



p is constant along a dipolar field line while q parameterizes the displacement parallel to the field: q = at 
the equator, q — > — oo as 9 — ► ix and q — ► +oo as 9 — ► 0. 

1 Inverse Transformation 

The inversion of Q — *- e - ; finding (r, 9) given (q,p) — involves the solution of a non-trivial equation. 
Substituting for r gives qp 2 = cos 9/ sin 4 6? while eliminating 9 leads to the polynomial expression 



Descartes's rule of signs states that ^ has exactly one positive, real root and, as has been previously noted 
[3], since J2J) is a quartic this root his an algebraic representation. To be useful in numerical models however 
such a solution has to be expressed in a computationally stable form. Define the auxiliary quantities 
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Then the positive real root of J2J is 
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Since ji > 1 this formulation eliminates the possibility of a catastrophic cancellation between terms. To 
complete the inversion recall that sin 2 9 = r/p. 

Although (J5J makes it possible, in principle, to express formulas in terms of either r and 9 or q and p, it 
is usually simpler to use the former representation. 

2 Coordinate Derivatives 

In this section we ignore the <j) coordinate. The partial derivatives of the dipolar coordinates with respect to 
the spherical polar coordinates are 
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The Jacobian is then 
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where 6(9) = \f\ + 3 cos 2 9. The derivatives of the spherical polar coordinates with respect to the dipolar 
coordinates are 
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A few second derivatives are also occasionally useful 
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3 Unit Vectors 

The expressions for the dipolar and spherical polar unit vectors in the alternate coordinate system are 
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r = -— q+ — p 9=-— Q-— P (14) 

3.1 First Derivatives of Unit Vectors 
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3.2 Some Second Derivatives 
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4 Metric; Differential Line, Area, and Volume Elements 
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Since the dipole system is orthogonal the only non-zero components of the metric are the diagonal elements 
(i.e., gij = for i ^ j). The associated scale factors (hi, where h 2 = go) are 
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With these we can immediately write down the differential elements 
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The Christoffel symbols of the second kind, as defined by Arfken £Q, are given by the formula 
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where g %% = 1/gu- In matrix form they are 
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5 Vector Operations 
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The differential operators can be derived from the metric tensor. In what follows / is a scalar, A and B arc 
vectors, and T is a tensor. 



5.1 Gradient 
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which is equivalent to 
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5.4 Scalar Laplacian 

V 2 / = 

which can also be written as 
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5.5 Vector Laplacian 
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5.6 Directional Derivative 
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5.7 Divergence of a Tensor 



V-T = q 

+P 

+4> 



3 sin 6 cos $ 3 

V • (T TC q + T pq p + T 09 0) + -^-(1 + cos 2 9)J qp + + cos 2 9)J PP + - 



rS 3 



„ ._ „ „ 7. 3sin0... 2 /!\-r- 6cos0._, o^-r- 
V • (T qp q + Tppp + T^pt/i) rs-U + cos 6)T qq (1 + cos 0)T P 



r6 3 



1 



r<5 sin 



— (1-3 cos" 0)T« 



-r - -r 3COS0 

v • (Tq^q + Tp^p + T000) — T+,0 + 



1 



rS sin 



r5 3 



(1-3 cos 2 0)T 0? 



Acknowledgments 

This work was supported by the Office of Naval Research. 



References 

[1] Arfkcn, G., Mathematical Method for Physicists, 3rd ed., San Diego, CA: Academic Press, pp. 158-163, 
1985. 

[2] Fatkullin, M. N., and Yu. S. Sitnov, Dipole coordinate system and some of its characteristics, Geomag- 
netism and Aeronomy, 12, 293-295, 1972, translated from Russian. 

[3] Kageyama, A., T. Sugiyama, K. Watanabe, and T. Sato, A note on the dipole coordinates, Computers 
and Geosciences, 32, 265-269, 2006. 

[4] Orens, J. H., T. R. Young, Jr., E. S. Oran, and T. P. Coffey, Vector operations in a dipole coordinate 
system, NRL Memo. Rep. 3984, Naval Research Laboratory, Washington, D.C., 1979. 



6 



